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Abstract 

We report on first measurements of low-mass electron-positron pairs in Pb-Au collisions at 
the SPS beam energy of 40 AGeV. The observed pair yield integrated over the range of invariant 
masses 0.2<m<l GeV/c'^ is enhanced over the expectation from neutral meson decays by a factor of 
5.9ibl.5(stat.)±1.2(syst. data)±1.8(syst. meson decays), somewhat larger than previously observed 
at the higher energy of 158 AGeV. The results are discussed with reference to model calculations 
based on tt'^'tt"— ^e+e~ annihilation with a modified /9-propagator. They may be linked to chiral 
symmetry restoration and support the notion that the in-medium modifications of the p are more 
driven by baryon density than by temperature. 
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reactions like p-Be, p-Au at 450 GeV, 
AGeV Q and Pb-Au at 158 AGeV 



According to finite temperature lattice QCD, strongly interacting matter will, at suf- 
ficiently high energy densities, undergo a phase transition from a state of hadronic con- 
stituents to quark matter, a plasma of deconfined quarks and gluons. At the same time, 
chiral symmetry, spontaneously broken in the hadronic world, will be restored. High-energy 
nucleus-nucleus collisions provide the only way to investigate this issue in the laboratory. 
Among the different observables used for diagnostics of the hot and dense fireball formed in 
these collisions, lepton pairs are particularly attractive. In contrast to hadrons, they directly 
probe the entire evolution of the fireball; the prompt emission after creation and the absence 
of any final state interaction conserve the primary information within the limits imposed by 
the space-time folding over the emission period. 

The CERES/NA45 experiment at the CERN SPS is the only dielectron spectrometer 
in the field of ultrarelativistic nuclear collisions, focused on the measurement of electron- 
positron pairs in the invariant mass range mee<l GeV/c^. In a series of systematic mea- 
surements, CERES has previously foimc^ronounced differences between proton-induced 

and nucleus-nucleus collisions like S-Au at 200 
^. While the superposition of known electromag- 
netic decays of produced neutral mesons can successfully account for the measured e^e~ 
mass spectra in the proton case, a strong excess of e~^e~ pairs above the expectation from 
meson decays has consistently been observed in the nuclear case, amounting for 158 AGeV 
Pb-Au collisions to a factor of 2.4±0.2(stat.)±0.6(syst. data)±0.7(syst. meson decays) in 
the mass interval 0.2<mee<l GeV/c^ (combined 1995/96 data 0,0,0]). The observation 
of this excess has initiated an enormous theoretical activity (see ja, LZf for recent reviews). 
There is a general consensus that one observes direct thermal radiation from the fireball, 
dominated by two-pion annihilation 7i~^TT~^p—>-e~^e~ with an intermediate p vector meson. 
The p is of particular relevance (more so than the other light vector mesons u and 0), due to 
its short lifetime of 1.3 fm/c and its direct link to chiral sjTiimetry restoration Q. Indeed, 
a proper description of the observed dilepton excess does require a strong modification of 
the intermediate p in the hot and dense medium. The two main theoretical alternatives 
for this modification are (i) 'Brown-Rho scaling^ith an explicit connection to the medium 
dependence of the chiral condensate , reducing the mass of the p below the 

vacuum value, and (ii) calculations of the p spectral density on the basis of p-hadron in- 
teractions ^ 12, spreading the width of the p above the vacuum value. Theoretically, 
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the modifications are more sensitive to the baryon density than to the temperature of the 
fireball. 

To critically examine the relative importance of baryon density and temperature, an 
experimental variation of these parameters seemed mandatory. In this Letter, we present 
first results [14] on the production of e^e^ pairs for Pb-Au collisions at the SPS beam energy 
of 40 AGeV, where a lower temperature and a higher baryon density (by a factor of 1.5) are 
expected as compared to 158 AGeV We indeed observe a strong excess of pairs again; 
the enhancement factor relative to the expectation from meson decays is now found to be 
5.9±1.5(stat.)±1.2(syst. data)±1.8(syst. decays), i.e. somewhat greater than at the higher 
beam energy. 

The CERES/NA45 experiment is optimized for the measurement of e+e^ pairs in the 
invariant mass range of tens of MeV/c^ up to 1 GeV/c^. The acceptance, azimuthally 
symmetric, covers the pseudorapidity region 2.11<riL<2.64: in the laboratory frame, i.e. a 
window of A?2=0.53 around mid-rapidity. A detailed description of the experiment can 
be found in jlfil . IitI 1^|. Here we summarize the most relevant features. Two Si drift 
detectors (SIDCl/2) 17j, located 10/14 cm downstream of a segmented Au target, provide 
an angle measurement of the charged particles and reconstruction of the vertex; at the same 
time, they serve as a powerful handle to recognize photon conversion- and Dalitz electron- 
positron pairs with their small opening angles. The basic discrimination between the rare 
electrons and the abundant hadrons is done with two Ring Imaging Cherenkov detectors 
(RlCHl/2) Q- For the run periods 1999/2000, the spectrometer was upgraded by 
the addition of a cylindrical Time Projection Chamber (TPC) [ISJ with a radial electric 
drift field, operating inside a new magnetic field configuration. As a consequence, the mass 
resolution is potentially improved from the 1995/96 value of 5.5% to about 2% in the region 
of the p/uj. Since the original magnetic field between the RICH's is no longer required, it is 
turned off, and the two detectors are now used as an integral unit, resulting in an improved 
electron efficiency (~0.9) and rejection power. 

The results reported in this Letter were obtained from the analysis of data taken during 
the SPS run period in the fall of 1999, the first on e~*"e~-pair production after completion 
of the upgrade. The centrality threshold, based on a lower limit of the integrated signal 
of the SIDCl detector (proportional to the total number of hits), was set to correspond to 
the 30% most central fraction of the geometrical cross section (same as for the 1996 data). 
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Since the new TPC readout was still being commissioned, only about 8 million events were 
taken. Related problems also limited the track efficiency of the TPC to about 0.43 and the 
mass resolution actually achieved to about 6% in the region of the p/u (inferred from an 
independent analysis of A and Kq production). 

The physics signal to be analyzed, e^e~ pairs with masses mee>0.2 GeV/c^, has an abun- 
dance of only 10"^ both relative to hadrons and to photons. The experimental challenges 
to minimize the associated background are therefore extremely high. The electron track re- 
construction and hadron suppression proceed in the following way. First the event vertex is 
reconstructed from 100 or more particle trajectories traversing both SIDC's. Track segments 
are then formed in the SIDC's. Independently, Cherenkov ring candidates are reconstructed 
according to a pattern recognition algorithm based on a two-step Hough-transformation for 
candidate search and a subsequent ring ffiting procedure Track segments from the 

SIDC's and RICH's are then matched, defining electron track candidates, and these can- 
didates are linked to the TPC and matched there to TPC track segments, independently 
preselected by dE/dx information characteristic for relativistic electrons. With a RICH ra- 
diator threshold of 7t/i=32, more than 99% of all charged hadrons are already rejected on the 
SIDC-RICH level. The remaining small pion contamination, due to accidental matches with 
fake-ring electron track candidates and to true matches with high-pt pions with asymptotic 
ring radii, is completely removed by the (momentum- dependent) dE/dx-cut in the TPC, 
a considerable improvement compared to the 1995/1996 analysis. Finally, the remaining 
electron and positron tracks are combined into pairs, and the type of pair is defined. 

The photon-related background is much more severe. Even with minimal material in the 
experiment, photons converting in the target and in SIDCl together with vr^-Dalitz decays 
still exceed the expected high-mass signal by a factor of about 10^, thus constituting the vast 
majority of the electron tracks. Although the characteristics are different (small pair opening 
angles and very low masses mee<0.2 GeV/c^), the limited track reconstruction efficiency and 
acceptance lead to a combinatorial high-mass background for events in which two or more 
of these low-mass pairs are only partially reconstructed. This is the central problem of 
the experiment. It is dealt with in the following way 14j. Since the inclusive electron 
Pt-spectra from low-mass pairs are considerably softer than those of pairs with mee>0.2 
GeV/c^, a reduction of the combinatorial background by more than a factor of 10 can be 
obtained by pairing only electron tracks with transverse momenta pt>0.2 GeV/c. Further, 



5 



a global selection on opening angles 0ee>35 mrad is applied to all pairs with complete 
tracks. Next, conversion- and Dalitz pairs with opening angles Bee<10 mrad which are not 
recognized as separate tracks with two individual electron rings in the RICH detectors are 
rejected by a correlated donhle-dE/dx cut in each of the two SIDC's 3|- Electron tracks 
which have a second SIDC-RICH electron candidate track within 70 mrad are also rejected. 
Finally, identified Dalitz pairs (mee<0.2 GeV/c^) are excluded from further combinatorics. 
Altogether, a rejection factor of more than 100 is obtained. 

The remaining combinatorial background is determined by the number of like-sign pairs. 
The net physics signal S is then obtained by subtracting the like-sign contribution from the 

total unlike-sign pairs as = A^+_ — 2(A^+_|_ ■ )^/^. The final sample after all cuts consists 

of 532 unlike-sign, 283 like-sign and 249±29 net unlike-sign pairs (signal-to-background ratio 
of 1/1) for masses mee<0.2 GeV/c^, and 1236 unlike-sign, 1051 like-sign and 185±48 net 
unlike-sign pairs (signal-to-background ratio of 1/6) for masses 0.2<mee<1.15 GeV/c^. 

The raw net yields require an absolute normalization and an absolute correction for re- 
construction efficiency for further discussion; an acceptance correction is not done. The 
normalization refers to the total number of events and the average number (Nch) of charged 
particles per event within the e^e^-pair rapidity acceptance. The charged particle rapidity 
density distribution is extracted from an independent analysis of the SIDC track segments; 
for the present centrality selection, one gets (A'"ch) = ((iA'"c/i/(i?7)-A?7=216-0. 53=115. The pair 
reconstruction efficiency e is determined by a Monte Carlo procedure in which simulated 
lepton pairs from neutral meson decays with known characteristics (see below) are embed- 
ded into real data and are then analyzed using the same software chain as used for the data; 
this 'overlay-MC provides the highest degree of realism, in particular for the description 
of event background and instrumental deficiencies. A value of (l/e)~^=0.036 is obtained 
as the weighted average over the (iAch/c??7-distribution (the only essential functional depen- 
dence), using as weights the uncorrected differential pair yield for mee<0.2GeV/c^. The 
primary reason for this low value (compared to 0.11 in 1996) is the 1999 status of the TPC 
readout, while the influence of the low-mass pair rejection cuts on the efficiency is found 
to be minor. This is conflrmed by the stability of the two net pair samples with masses < 
and > 0.2GeV/c^ as a function of the rejection steps over a wide dynamical range of the 
signal-to-background ratio jl^ . 

The normalized and efficiency-corrected mass spectra for the total unlike-sign pairs and 
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the like-sign pairs are shown in Fig. The statistical bin-to-bin fluctuations of the mass- 
differential net signal can be reduced by about a factor of v^, if the shape of the like-sign 
combinatorial background is determined with much better accuracy. Somewhat different 
from previous CERES analyses, this is done here by using the like-sign combinatorial back- 
ground after the first two rejection steps (only pt- and global Qee- cut) which contains 
about 10 times more entries than the final background. This 'high-statistics' background is 
included in Fig. ^ as a histogram, rescaled as to give the same integral as the final back- 
ground; the two are consistent within errors. The pair transverse momentum spectra are 
treated similarly. 

The resulting net invariant- mass spectrum of e^e^ pairs is shown in Fig. |21 The errors 
attached to the data points are purely statistical. The systematic errors contain contri- 
butions from the normalization procedure and the efficiency correction, amounting to an 
overall error of 16% for masses mee<0.2 GeV/c^. In the high-mass region mee>0.2 GeV/c^, 
additional errors of about 20% arise from remaining ambiguities in the rejection cuts and 
from the statistical error of the background rescaling, raising here the overall error to 26%. 
All values are quadratically small compared to the individual statistical errors of the data 
points and are therefore not drawn. Integration of the mass spectrum results in total nor- 
malized e"'"e~-pair yields per event per charged particle of [7.6±0.8(stat.)±1.2(syst.)] -10"^ 
for mee<0.2 GeV/c^, and [5.6±1.4(stat.)±l.l(syst.)] -lO^^ for mee>0.2 GeV/c^. The back- 
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FIG. 1: Mass spectra of unlike-sign and like-sign pairs (circles). A high-statistics like-sign spectrum 
is shown as a histogram (see text). The peak at mee>0.4 GeV/c^ is due to the single-electron cut 
Pt>^.2 GeV/c. 
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ground rescaling used to improve the differential spectrum has, of course, no relevance for 
the integrals. Their statistical errors are those of the primary samples quoted above, and 
the systematic error of the high-mass integral mee>0.2 GeV/c^ is only 20% (instead of 26%). 
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FIG. 2: Inclusive e~^e~ mass spectrum, compared to the hadron decay cocktail (thin solid; individ- 
ual contributions thin dotted) and to theoretical model calculations 19] based on vr+vr^ annihilation 
with an unmodified p (thick dashed), an in-medium dropping p mass (thick dashed-dotted) and an 
in-medium spread p width (thick solid). The model calculations contain the cocktail, but without 
the p to avoid double counting. The low-mass tail of the cocktail p is due to the inclusion of a 7r"^7r~ 
phase space correction. The (weaker) tails of the to and (p are caused by electron bremsstrahlung. 

Fig. 121 also contains a comparison to the electromagnetic decays of the produced neutral 
mesons. The evaluation of this 'hadron decay cocktail' follows our previous procedure, 
described in detail for proton-induced reactions in (211 and extended to the Pb-Au case 

nn 

in 1^ y| , including the empirical systematics of hadron yields (statistical model values where 
unmeasured) and of mj-distributions (influenced by radial flow). The simulations are subject 
to the same cuts as the data; resolution broadening as inferred from hadron data (see above) 
and the effects of bremsstrahlung from the electrons traversing a total detector material of 
5.8% of a radiation length before the TPC are included. For masses mee<0.2 GeV/c^, 
the spectrum is dominated by 7r°- and rj Dalitz decays, and good agreement with a ratio 
data/decays of 1.18±0.13(stat.)±0.19(syst. data)±0.09(syst. decays) is found. For masses 
mee>0.2 GeV/c^, however, a structureless continuum much above the hadron-decay level is 
seen, extending all the way up to the region of the 0. The enhancement factor data/decays 
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is 5.9±1.5(stat.)±1.2(syst. data)±1.8(syst. decays) in this region, i.e. larger than the value 
at the higher beam energy of 158 AGeV with a statistical significance of 1.8a (data errors 
added in quadrature); the systematic errors 'decays' of the decay cocktail, dominated by 
uncertainties in the branching ratios and transition form factors, essentially drop out in the 
comparison. 

The invariant differential spectra in pair transverse momentum p^'^ are plotted in Fig. El 
separately for the two mass regions mee< 0.2 and >0.2 GeV/c^. The steep cut-off in the dis- 
tribution for pI^kOA GeV/c (left part) reflects the cut of the single-electron pt>0.2 GeV/c. 
Comparing data and hadron decays, good agreement is seen again for masses mee<0.2 
GeV/c^, while the strong excess above the decay expectation reappears for the high-mass 
region, mostly localized at low pf<Q.Q GeV/c as observed before at 158 AGeV 0,0,01. 

To illustrate most directly the relevance of the new results. Figs. El and El also contain 
the predictions from theoretical model calculations Q|, based on tt+tt^ annihilation as 
discussed in the introduction and including the hadron decay contribution. The /?-pro pag ator 
is treated in three ways - vacuum p, modiflcations through 'Brown-Rho scaling' ,yiiil, 
and modiflcations through p-hadron interactions ^ Q, Q|. With reference to the mass 
spectrum in Fig. |2l the data clearly rule out an unmodifled p, but agree with the two in- 
medium scenarios within errors (without being able to distinguish between them). The 
pair transverse momentum spectra in Fig. El are also reasonably well described, but the 
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FIG. 3: Invariant pair-pf*^ spectra for masses mee<0.2 GeV/c^ (left) and mee>0.2 GeV/c^ (right), 
compared to the hadron decay cocktail and to theoretical model calculations Q|. The cuts, labels 
and line codes are identical to those in Fig. [21 
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peculiar rise (in contrast to the shape of the hadron-decay spectra) appears here as 

a basic feature of the pion-annihilation process without much discrimination power towards 
in-medium effects. The integral yield of the model calculations (thick solid line in Figs. 2 
and 3) corresponds to a 'theoretical' enhancement factor of 4.2 at 40 AGeV, compared to 
2.4 at the higher beam energy; the difference is consistent with the observed trend. 

We draw the following conclusions. The e"'"e~-pair yield for m>0.2GeV/c^ observed at 40 
AGeV is enhanced over the expectation from neutral meson decays. Compared to 158 AGeV, 
the enhancement factor may even be larger, with a significance of 1.8cr. However, even if it 
would be the same, that would still be remarkable vis-a-vis the changes in conditions, i.e. 
less multiplicity and less temperature at the lower energy. The only quantity rising is the 
baryon density. Therefore, on very general grounds and independent of the details of any 
model calculations, the sole existence of a strong enhancement at 40 AGeV is consistent with 
the theoretical description of vr+vr" annihilation and a p propagator dominantly influenced 

],Q,y,Qiii[y,QQ- 



by baryon density rather than temperature 



It is also consistent 



with theoretical results that contributions from qq annihilation during the QGP phase are 
minor 
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